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Table  1.  Carbohydrate  Information  Density.  The  information 
density  of  biopolymers  can  be  considered  the  number  of 
different  molecules  (N)  that  can  be  made  from  an  oligomer  of  a 
given  size  (SJ.  For  nucleotides,  W  =  4®  and  for  peptides  N  =  20®, 
but  for  carbohydrates  the  equation  is  far  more  complex. 

Number  of  Different  Oligomers 


Oligomer 

Size 

Nucleotides 

Peptides 

Carbohydrates 

1 

4 

20 

20 

2 

16 

400 

1,360 

3 

64 

8,000 

126,080 

4 

256 

160,000 

13,495,040 

5 

1,024 

3,200,000 

1,569,745,920 

6 

4,096 

64,000,000 

192,780,943,360 

A.  Introduction  and  Objectives 

Self-assembly  provides  an 
attractive  route  from  simple  building 
blocks  to  hierarchical  nanomaterials 
whose  sophistication  and  function  rival 
those  found  in  living  systems.  Nature 
uses  the  cooperative  and  multivalent 
interactions  between  carbohydrates 
and  the  proteins  that  bind  them,  termed 
lectins,  to  direct  the  formation  of  larger 
aggregates,  such  as  tissues,  organs, 
and  ultimately  organisms.  Despite  their 
increased  information  density  relative 
to  other  recognition  elements  (DNA 
base  pairs,  peptides,  etc.)  currently 

used  to  assemble  hierarchical  _ 

nanosystems  (Table  1),  carbohydrates 

have  not  been  adopted  by  researchers  for  organizing  synthetic  assemblies  because  natural 
lectins  are  too  complex  to  use  in  the  context  of  nanotechnology,  and  the  intricacies  of 
oligosaccharide  structure  makes  designing  specific  hosts  that  can  discriminate  between 
different  monosaccharides  particularly  difficult.  While  there  is  significant  interest  in  developing 
synthetic  lectins,  this  task  remains  challenging  because  monosaccharides  contain  several 
hydroxyl  functions  and  multiple  stereocenters  per  residue,  and  they  typically  occur  as  a 
constantly  fluctuating  mixture  of  5-  and  6-membered  rings.  However,  the  multivalency  and 
cooperativity  inherent  to  saccharide  binding  could  add  significant  capabilities  to  synthetic 
nanotechnology.  Both  carbohydrates  and  lectins  typically  exist  as  multimeric  assemblies,  and 
the  high  affinity  that  arises  from  multivalent  interactions  leads  to  the  error-correcting  abilities  that 
are  characteristic  of  biological  self-assembly.  Moreover,  binding  in  lectins  can  change  as  a 
function  of  the  cooperativity  that  arises  from  the  interplay  of  two  or  more  interactions  with  the 
same  receptor,  so  that  the  system  as  a  whole  behaves  differently  than  the  mere  multiplication  of 
individual  interactions  would  predict.  The  goal  of  the  research  program  was  to  utilize 
oligosaccharide-synthetic  lectin  binding  pairs  to  direct  the  assembly  of  self-organized, 
hierarchical  materials  and  to  employ  the  binding  properties  of  these  structures  to  create 
synthetic  nanosystems  for  information  storage  and  processing.  To  that  end,  we  pursued  three 
synergistic  objectives:  (1)  develop  saccharide-  synthetic  lectin  binding  pairs  that  bind  specifically 
and  strongly  through  noncovalent  interactions,  (2)  utilize  the  information  encoded  by  the 
oligosaccharide  structure  to  design  advanced  hierarchical  materials  for  molecular  logic 
operations,  (3)  employ  cooperative  and  multivalent  binding  interactions  between  carbohydrates 
and  synthetic  lectins  to  drive  mesoscopic  assembly  in  solution,  in  polymers,  and  on  surfaces. 
The  major  questions  we  asked  were:  How  can  we  develop  synthetic  lectins  that  discriminate 
between  monosaccharides  that  differ  by  only  the  orientation  of  a  single  hydroxyl?  How  can  the 
conformational  and  binding  dynamics  inherent  to  saccharides  be  used  in  molecular  logic?  Can 
synthetic  lectins  employ  multivalency  and  cooperativity  to  control  the  structure  and  dynamics  of 
nanomaterials  in  three  dimensions? 
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B.  Progress  Towards  Major  Research  Goals 

During  the  award  period,  the  majority  of  the  project  goals  were  accomplished,  and 
significant  progress  was  made  towards  all  research  goals.  We  have  developed  a  new, 
mannose-specific  synthetic  lectin  scaffold  that  achieves  specificity  through  cooperative  and 
multivalent  pathways  (Section  C.1.).  The  ability  of  this  molecule  to  switch  between  targets  and 
the  complex  molecular  logic  that  drives  this  switching  has  been  determined,  and  this  is  a 
successful  demonstration  of  one  of  the  major  project  goals:  demonstration  of  a  synthetic, 
carbohydrate  based  information  processing  system.  Another  major  aim  of  the  proposed  work 
was  to  optimize  new  surface  chemistries  for  immobilizing  carbohydrates.  To  this  end,  we  have 
developed  the  first  chemistry  to  covalently  pattern  graphene  using  force-accelerated  Diels-Alder 
reactions  (Section  C.2.),  and  this  chemistry  is  compatible  with  carbohydrates.  Subsequently, 
we  used  massively-parallel  tip  arrays  to  deliver  light  sight-specifically  to  a  surface  and  thereby 
initiate  photochemical  brush  polymerizations.  This  reaction  methodology  was  used  to  create 
carbohydrate-containing  brush  polymers  with  binding  affinity  that  is  orders-of-magnitude  greater 
than  conventional  glycan  microarrays  (Section  C.3.).  Because  of  the  critical  role  that  force  plays 
in  surface  reactions,  we  have  systematically  and  quantitatively  investigated  how  force 
accelerates  reactions  on  surfaces,  and  developed  a  new  predictive  model  explaining  previously 
contentious  observations  related  to  the  force-accelerated  Huisgen  reaction  (Section  C.4.), 
which  we  have  used  previously  to  immobilize  carbohydrates.  Finally,  we  have  continued  to 
explore  the  dynamics  of  activated  suspensions  consisting  of  Au-Pt  nanorods,  and  we  have 
developed  a  model  to  describe  their  movement  around  microspheres  (Section  C.5.),  and  these 
studies  will  move  us  towards  the  goal  of  creating  integrated  carbohydrate-containing 
nanosystems  with  complexity  of  natural  information  carrying  networks. 
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C.  Results  and  Discussion 

C.1.  Synthetic  Lectin  Design.  Natural  saccharide-binding  proteins,  including  lectins  and 
periplasmic  substrate-binding  proteins,  use 
water  desolvation,  hydrogen  bonding  (H- 
bonding),  and  C-H***tt  interactions  to 
selectively  recognize  glycans  that  may  differ 
only  by  the  orientation  of  a  single  hydroxyl 
group  to  achieve  binding  affinities,  KgS,  as 
high  at  10®  M\  Selective  carbohydrate 
recognition  with  artificial  receptors  remains  a 
major  area  of  investigation  because  of  the 
challenge  of  differentiating  between 
molecules  with  subtle  structural  differences, 
their  ability  to  reveal  fundamental  aspects  of 
saccharide  binding,  and  their  potential 
applications  in  disease  detection,  therapy,  or 
catalysis.  Mannose  is  a  particularly 
interesting  monosaccharide  target  because  it  is  involved  in  signaling  for  several  cancers,  and  as 
a  consequence  developing  mannose  specific  synthetic  receptors  could  reveal  how 
carbohydrates  are  employed  for  communication  purposes.  The  challenge  with  mannose  as  a 
target  for  any  of  these  applications  is  that  it  is  an  epimer  of  glucose  -  the  most  prevalent 
monosaccharide  in  biology  -  and  is  thus  very  difficult  to  design  receptors  that  can  distinguish 
between  molecules  that  differ  only  by  the  orientation  around  a  single  stereocenter.  However, 
synthetic  carbohydrate  receptors  with  increased  binding  affinity,  expanded  substrate  scope 
beyond  all-equatorial  glycosides,  and  increased  selectivity  continue  to  be  a  major  research  goal 
because  of  the  new  capabilities  they  could  provide  in  terms  of  biological  activity,  sensing,  and 
hierarchical  assembly. 
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Fig.1  Structures  of  1 
whose  binding  with  1 
numbering  scheme. 
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Fig.  2  The  equilibria  present  in  a  chloroform  mixture  of  1  and  p-Man  at  25°C.  Upon  increasing  the  p-Man:1  ratio,  K2 
determines  the  dominant  species  in  solution.  If  instead,  p-Man:1  decreases,  K3  controls  the  equilibrium  mixture. 


To  address  these  challenges,  we  pursued  an  alternate  design  approach  towards 
synthetic  saccharide  receptors  that  involves  preparing  a  flexible  host  that  does  not  possess  rigid 
preorganization  and  as  a  result  is  capable  of  binding  carbohydrates  through  pathways  that  arise 
from  conformational  rearrangements  and  positive  homotropic  cooperativity.  Homotropic 
cooperativity  -  whereby  an  initial  association  of  a  target  substrate  induces  conformational 
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restrictions  that  enhance  further  binding  of  the  same  substrate  -  has  been  employed  previously 
in  synthetic  receptors  to  increase  binding  strength  towards  targets  such  as  diacids  and  syn- 
diols.  However,  synthetic  receptors  that  utilize  homotropic  cooperativity  remain  rare,  and  few 
examples  exist  of  hosts  that  rely  upon  cooperativity  to  enhance  saccharide  binding,  despite  the 
fact  that  cooperativity  and  multivalency  are  ubiquitous  elements  of  carbohydrate  recognition  in 
biology.  The  advantages  that  could  arise  with  carbohydrate  receptors  that  employ  positive 
allosteric  cooperativity  include  (1)  shorter  synthetic  sequences  because  specificity  and  affinity 

are  the  direct  result  of  allostery,  (2)  increased 
solubility  associated  with  receptor  flexibility,  and 
(3)  employing  the  same  binding  mechanisms  as 
natural  saccharide  receptors  could  provide 
fundamental  insights  into  the  complex 
carbohydrate  recognition  motifs  that  are  prevalent 
in  nature.  To  this  end,  we  have  designed  receptor 
1  that  contains  aminopyrrolitic  groups  appended 
to  a  flexible  scaffold  with  eight  freely  rotating 
methylene  groups  and  a  freely  rotating  biphenyl 
bond  (Figure  1).  Thus,  1  was  designed  to 
dynamically  explore  thermodynamic  and 
conformational  space  and  confirm  that  increased 
receptor  flexibility  can  induce  specificity  for 
carbohydrate  guests  through  allostery  despite  a 
higher  entropic  penalty  experienced  in  an  initial 
association  step.  Through  the  use  of  NMR 
titrations,  we  found  that  1  achieves  excellent 
selectivity  for  octyl  mannosides  through  two 
allosteric  cooperative  pathways  (Figure  2)  with  an 
overall  selectivity  as  high  as  16.8:1  a-Man:a-Gal 
in  chloroform.  Moreover,  the  selectivity  is  directly 
dependent  on  pyranoside  concentration,  where  1 
binds  preferentially  to  p-GIc  at  low  concentration 
(<0.3  mM),  then  binds  a-  and  p-Man  at  higher 
concentrations.  To  the  best  of  our  knowledge,  1  is 
the  first  synthetic  carbohydrate  receptor  that  (1) 
relies  on  cooperativity  to  increase  selectivity  and 
(2)  whose  selectivity  switches  with  saccharide 
concentration. 

Pyranoside  binding 

Receptor  1  was  prepared  in  a  five  step  synthetic 
sequence  in  a  64%  overall  yield,  and  the  binding 
affinity  and  selectivity  of  1  for  eight  octyl 
pyranosides  (Figure  1),  which  were  chosen  as 
guests  because  they  are  common  terminal 
residues  found  on  cell  surface  glycoconjugates 
and  are  used  as  standards  for  assaying  synthetic 
receptor  selectivity  were  investigated  by  variable  temperature  NMR  titrations  in  CDCI3. 
Chloroform  facilitates  the  binding  studies  of  new  carbohydrate  receptors  because  the  solvent 
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Fig.  3  a)  NMR  (600  MHz,  CDCI3,  25^) 
spectra  obtained  after  the  incremental  addition  of 
a  10.4  mM  solution  of  p-Man  to  a  58.6  pM 
solution  of  1,  with  dashed  lines  illustrating  the 
induced  changes  in  chemical,  shifts,  b)  The 
chemical  shift  dependences  of  of  1  (58.6  pM)  is 
plotted  as  a  function  of  molar  equivalents  of  each 
of  the  eight  added  pyranosides  from  the  NMR 
titrations,  c)  The  fittings  of  the  chemical  shift 
changes  in  response  to  the  addition  of  p-Man  to  a 
model  containing  K^  (dashed  line)  and  /Ci  +  K2 
(solid  line). 


does  not  compete  for  H-bonding  between  host  and  guest,  thereby  enhancing  polar  noncovalent 
bonds.  These  titrations  revealed  that  the  selectivity  of  1  for  octyl  mannosides  arises  through  a 
cooperative,  multistep  equilibrium  (Figure  2).  In  solution,  1  exists  as  a  dimer,  I2,  that  comes 
apart  to  form  1:1  complexes  with  the  eight  pyranosides  with  little  selectivity.  Upon  altering  the 
lipyranoside  ratio  from  1:1,  two  new  mannoside-selective  binding  modes  emerge  that  are  both 
the  result  of  positive  cooperativity  at  25°C.  When  the  concentration  of  mannoside  is  increased, 
1:2  receptonpyranoside  complexes  form  only  with  a-Man  (1:a-Man2)  and  3-Man  (1:3-Man2).  If 
instead,  the  concentration  of  1  increases  with  respect  to  pyranoside,  a  2:1 
receptor:monosaccharide  complex  forms  only  with  p-Man  (l2:P-Man).  To  understand  the  origin 
of  the  preferential  binding  of  mannosides  by  1 ,  these  binding  stoichiometries  and  the  structures 
of  12:3-Man,  1:a-Man2,  and  1:3-Man2  were  established  through  a  variety  of  ID  and  2D  NMR 
methods,  and  all  values  of  Ku  K2,  K3,  AhP  and  AS°  were  determined  for  each  of  these  equilibria. 


1:2  Receptonpyranoside  binding 

The  binding  of  pyranosides  was  first  studied  under  conditions  where  receptor  1  was  maintained 
at  a  low  concentration  (<70  pM)  to  minimize  the  contribution  of  the  K^imer  and  K3  equilibria 
(Figure  2).  Once  the  values  of  K^  and  K2  were  determined,  their  values  were  held  invariant  in 
subsequent  peak  shift  fittings,  which  facilitated  the  quantification  of  the  other  KgS.  The 
incremental  addition  of  a  10-12  mM  solution  of  each  of  the  eight  octyl  pyranoside  to  a  dilute 
solution  of  1  (58.6  pM)  induced  changes  in  the  ^Fl  NMR  chemical  shifts  corresponding  to  the 
protons  of  1,  owing  to  a  fast  exchange  between  bound  and  unbound  substrates  on  the  NMR 
timescale  (Figure  3A).  Notably,  the  change  of  mannoside  proton  chemical  shifts,  A6,  was 
significantly  greater  than  the  other  pyranosides,  suggesting  that  the  overall  binding  free  energy, 
AG°,  was  greater  for  mannosides  than  other  pyranosides  (Figure  3B).  When  the  observed 
chemical  shifts  were  plotted  against  the  molar  equivalents  of  pyranoside  (Figure  3C),  the 
observed  perturbations  for  the  mannoside  protons  follow  a  sigmoidal  pattern,  indicating  that 
multiple  equilibria  in  addition  to  1:1  complexation  are  occuring  in  the  titration.  The  titration 
induced  peak  shifts  occurring  upon  addition  of  both  anomers  of  glucose  (a/p-GIc),  galactose 
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(a/(3-Gal),  and  N-acetylglucosamine  (a/(3-GlcNAc)  produced  no  such  sigmoidal  curve  and 
instead  follow  the  hyperbolic  shape  of  a  1:1  binding  isotherm  (Figure  3B).  Since  the 
concentration  of  pyranosides  was  intentionally  kept  low,  the  saturation  region  of  the  binding 
isotherm  was  only  reached  for  the  mannosides  (Figure  3C),  thus  K^  could  not  be  determined  for 
the  other  pyranosides  from  these  experiments,  although  we  were  able  to  obtain  these  values  by 
titrating  1  into  solutions  of  pyranoside. 

The  chemical  shift  changes  that  occur  because  of  the  interaction  of  1  with  a-  and  3-Man 
were  subjected  to  a  global  nonlinear  fitting  analysis  with  a  model  incorporating  K'dimer,  K^,  and  K'a, 
and  satisfactory  fits  to  the  A6s  were  obtained  to  provide  macroscopic  KaS  indicating  the 
presence  of  1:2  receptor:mannoside  complexes  1:a-Man2  and  1:p-Man2.  For  an  allosteric 
receptor  possessing  two  equivalent  binding  sites,  the  experimentally  measured  KgS  must  be 
corrected  for  the  existence  of  two  identical  1 :1  intermediates  to  obtain  microscopic  association 
constants  that  accurately  describe  the  association  of  each  binding  site.  Since  there  are  two 
identical  pathways  by  which  a/p-Man  can  associate  with  1  to  form  a  1:1  complex,  the 
macroscopic  Kg  for  the  first  association  process  was  divided  by  2  to  obtain  microscopic 
association  constant  K^.  Likewise,  as  there  are  two  identical  pathways  for  dissociation  in 
1  :Man2,  the  macroscopic  Kg  for  the  second  association  process  was  multiplied  by  2  to  obtain 
microscopic  association  constant  K'2  (Table  2).  The  interaction  parameter,  a,  is  the  ratio  of  a 
microscopic  Kg  in  the  presence  of  cooperativity  over  its  value  in  the  absence  of  cooperativity,  or 
reference  Kg,  and  is  a  quantitative 


measure  of  cooperativity.  Since  both 
binding  sites  are  identical,  the  reference 
Kg  is  equivalent  to  and  thus  a  =  K'2  / 
K^.  An  a  value  greater  than  1  indicates 
that  the  first  association  event  enhances 
the  second,  i.e.  positive  cooperativity. 


Table  2.  Binding  constants  {KaS,  25 °C)  and  thermodynamic 
parameters  (AH°,  AS°)  associated  with  the  1 :1  (upper 
value),  1:2  (middle  value),  and  2:1  (lower  value) 
1:pyranoside  binding  stoichiometries  determined  by  NMR 
titrations  and  the  intrinsic  median  binding  concentration 
BC5o°  calculated  from  K'dimer,  Ku  Kz,  and  K3  for  each 
pyranoside  interacting  with  1  in  CDCI3  at  25 °C. _ 


Likewise,  the  binding  is  negatively 
cooperative  if  a  is  less  than  1,  and  the 
binding  is  non-cooperative  if  a  is  equal  to 
1.  From  an  analysis  of  the  microscopic 
K'gS,  a-Man  and  p-Man  were  both  found 
to  possess  a  values  of  13.7  and  7.6 
respectively  with  receptor  1 ,  indicating  a 
high  degree  of  positive  cooperativity 
occurs  between  1  and  both  mannosides 
(Table  3). 

The  thermodynamic  origin  of  this 
positive  allosteric  cooperativity  in  the 
formation  of  1:p-Man2  was  investigated 
by  determining  AH°  and  AS°  associated 
with  each  binding  step.  The  NMR 
titrations  between  1  and  p-Man  were 
repeated  at  20,  15,  and  10°C,  and  K) 
and  K'2  values  for  each  temperature  were 
obtained  and  subjected  to  Van’t  Hoff 
analyses  to  provide  the  thermodynamic 
parameters  associated  with  each  binding 
event.  Notably,  the  AH°s  of  association 


Glycoside 

Log  7^ 

AH° 

(kcal  mol'^) 

AS° 

(e.u.) 

BCso'' 

(mM) 

a-Glc 

2.75  +  0.04 

-12.5  +  0.3 

-29  +  1 

a 

- 

- 

1.86  +  0.16 

a 

- 

- 

P-Glc 

3.16  +  0.01 

-16.5  +  0.1 

-41  +  1 

a 

- 

- 

0.70  +  0.02 

039" 

-6.6  +  0.4 

-21  +  1 

a-Man 

2.57  +  0.19 

c 

c 

3.71+0.10 

c 

c 

0.46  +  0.05 

a 

- 

- 

p-Man 

2.46  +  0.31 

-20.5  +  0.8 

-57  +  3 

3.34  +  0.11 

-20.2+1.2 

-52  +  4 

0.73  +  0.10 

2.45  +  0.09 

-11.0  +  0.4 

-26  +  2 

a-Gal 

2.18  +  0.02 

-13.1+0.5 

-34  +  2 

a 

- 

- 

7.74  +  0.31 

a 

- 

- 

P-Gal 

2.59  +  0.03 

-15.4  +  0.8 

-40  +  3 

a 

- 

- 

2.74  +  0.18 

a 

- 

- 

a-GlcNAc 

2.53  +  0.02 

-11.4  +  0.2 

-26  +  1 

a 

- 

- 

3.18  +  0.14 

a 

- 

- 

P-GlcNAc 

2.65  +  0.05 

-11.8  +  0.5 

-27  +  4 

a 

- 

- 

2.37  +  0.26 

a 

- 

- 

“  No  detectable  binding.  **  Estimated  by  extrapolation  of  a  van’t  Hoff  plot 
with  the  assumption  ACp  =  0.  Not  determined. 
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for  the  first  and  second  binding  step  are  identical,  -20.5  ±  0.8  and  -20.2  ±  1.2  kcal  mol'^ 
respectively,  suggesting  that  both  equivalents  of  (3-Man  bind  1  with  an  identical  number  of 
noncovalent  interactions  and  that  the  mannosides  do  not  interact  with  each  other  when  bound  to 
1.  A  comparison  of  the  corresponding  AS°  values  indicates  a  decrease  in  the  unfavorable  AS° 
occurs  in  the  second  binding  step  compared  to  the  first,  -52  ±  4  and  -57  ±  3  e.u.  respectively. 

The  geometries  of  1:(3-Man  and  1:(3-Man2  were  determined  by  ID  and  2D  NMR  and 

computational  modeling  to  understand  the 
structural  origin  of  the  cooperative  binding. 
The  NMR  spectrum  of  a  solution  of  1  (1 
mM)  and  (3-Man  (2  mM)  in  CDCI3  exhibits 
averaged  signals  for  all  protons  at  25  °C 
because  of  a  fast  exchange  of  the  signals 
corresponding  to  1,  (3-Man,  1:(3-Man2,  and 
1:3-Man.  Upon  cooling,  the  NMR 

resonances  corresponding  to  1,  except  H'^, 
decoalesce  into  two  sets  of  peaks  at  -40  °C 
(Figure  4),  indicating  a  partial 
desymmetrization  occurring  with  receptor  1. 
No  decoalescence  was  observed  for  the 
resonances  corresponding  to  (3-Man  signals 
upon  cooling  the  solution  (Figure  4), 
suggesting  that  both  equivalents  of  (3-Man 
occupy  identical  chemical  environments  in  the 
1:P-Man2  complex.  In  the  absence  of  3-Man, 
the  signals  corresponding  to  1  in  the  NMR 
spectrum  do  not  decoalesce  at  low 
temperature,  but  rather  are  broadened.  Thus, 
the  new  signals  that  arise  upon  cooling  of  the 
mixed  sample  can  be  attributed  to  a  complex 
consisting  of  1  and  p-Man.  A  comparison  of 
the  ^Fl  NMR  spectra  at  -63 °C  with  varying 
ratios  of  1:3-Man  revealed  that  the  broad 
signals  corresponding  to  free  1  disappear 
when  more  than  two  molar  equivalents  of  p- 
Man  are  present  in  solution,  further  confirming 
the  1 :2  stoichiometry  of  the  complex.  The  two 
sets  of  signals  observed  in  the  NMR 
spectrum  at  -63  °C  (Figure  4)  for  the  receptor 
and  only  one  for  the  mannoside,  indicate  that 
both  pyranosides  are  bound  to  two 

aminopyrrolitic  arms  of  1  and  are 

symmetrically  equivalent,  resulting  in  a  1:P-Man2  complex  that  exhibits  C2  symmetry.  A 
configuration  where  each  mannose  binds  to  two  aminopyrrolitic  arms  on  the  same  aromatic  ring 
would  render  both  H'^  protons  of  receptor  1  symmetrically  equivalent.  With  NMR  providing  a 
general  understanding  of  how  p-Man  sits  within  1,  a  molecular  mechanics  (AMBER*)  Monte 
Carlo  conformational  search  was  utilized  to  model  the  binding  geometry  of  the  1:p-Man 
complex.  The  conformational  searches  yielded  only  one  minimum  energy  structure  for  1:(3-Man, 
which  was  further  optimized  using  density  functional  theory  (B3LYP/6-31G-i-(d)).  The  resulting 
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Fig.  5  Energy  minimized  structure  (DFT,  B3LYP/6- 
31+(d))  for  1:p-Man  from  an  initial  binding  geometry 
obtained  by  an  AMBER*  mixed  low-mode/torsional 
Monte  Carlo  conformation  search.  Intermolecular  H- 
bonds  are  denoted  by  black  dashed  lines  and  C-H  '  tt 
interactions  are  denoted  by  orange  dashed  lines. 
Three  H-bonds  occur  between  the  hydrowl  group  of 
C2  and  both  amino  H®  and  one  pyrrole  H^protons  of 
the  receptor.  Additionally,  a  C-H  tt  interaction 
between  the  phenyl  ring  and  and  two  more  H- 
bonds  between  the  hydroxyl  group  of  C3  and  an 
amino  H®  and  pyrrole  proton  of  the  receptor  were 
observed.  Side  views  b)  parallel  and  c)  perpendicular 
to  the  biphenyl  linkage  of  the  energy  minimized 
structure  (DFT,  B3LYP/6-31+(d))  for  C2  symmetrical 
1:p-Man2.  The  biphenyl  dihedral  is  denoted  by  (p. 


calculated  structure  (Figure  5A)  is  in  excellent  agreement  with  the  experimental  1D  NMR 
data. 

The  origin  of  the  cooperativity  in  receptor  1  was  examined  computationally  by  varying 
the  dihedral  about  the  biphenyl  bond  {cp,  Figure  5B/C)  as  this  is  the  only  dynamic  element  of  the 
receptor  that  is  shared  by  both  binding  sites.  While  desolvation  contributes  significantly  to  AS°  it 
is  not  the  source  of  cooperativity  as  both  binding  sites  and  mannosides  are  solvated  identically. 
A  comparison  of  the  DFT  (B3LYP/6-31G+(d))  minimized  structures  corresponding  to  1,  1:3- 
Man,  and  1:3-Man2  revealed  cp  values  of  40.2°,  38.5°,  and  39.1  °  respectively,  indicating  that  the 
contribution  of  torsional  strain  towards  the  observed  cooperativity  is  negligible  because  (p  in  the 
bound  state  is  close  to  the  preferred  cp  of  unbound  1 .  More  likely,  the  presence  of  a  bound  p- 
Man  restricts  the  rotation  about  cp  and  incurs  an  entropic  penalty  that  is  only  paid  in  the  first 
association. 

1:1  and  2:1  Receptor :pyranoside  binding 
The  Ki  and  K3  of  the  binding  of  1  to 
pyranosides  were  determined  by  titrating  a 
solution  of  1  (53.0  -  62.5  mM)  into  a  1.0  mM 
CDCI3  solution  of  each  octyl  glycoside  until  a 
large  excess  of  1  was  reached.  Addition  of  1 
produced  considerable  changes  in  the 
NMR  resonances  of  the  receptor  and  all  eight 
pyranosides  (Figure  6A).  For  each  non- 
mannoside,  was  determined  by  subjecting 
all  resolvable  resonances  from  each  titration 
at  25°C  to  a  global  nonlinear  fitting  analysis 
with  a  model  combining  and  Kdimer  (13.0  ± 

0.5  M'^)  (Table  3).  Likewise,  when  the 
observed  upfield  chemical  shift  changes 
associated  with  a-Man  were  fit  to  a  model 
incorporating  Ksmex  and  the  previously 
determined  values  of  Ki  and  K'2,  excellent  fits 
for  the  peak  shifts  were  obtained  (Table  3). 

However,  the  peak  shifts  for  the  3-Man 
titration  could  only  be  fit  accurately  when  the 
formation  of  a  complex  consisting  of  two 
molecules  of  1  and  one  molecule  of  p-Man 
(l2:P-Man,  Figure  2)  was  considered  (Figure 
6B).  The  observed  A5s  were  fit  to  a  model 
consisting  of  Ks^ex,  ^1,  K2,  and  K3,  with  every 
binding  constant  except  K3  held  invariant,  to 
determine  K3  for  p-Man.  The  values  of  log  K^ 
for  binding  between  1  and  all  pyranosides 
range  from  2.5  to  3.3  with  the  highest  values 
associated  with  a-  and  p-Man  and  3-Glc.  Importantly,  little  difference  exists  between  the  KiS  for 
the  eight  pyranosides,  indicating  that  little  selectivity  occurs  in  the  first  binding  event.  In  contrast 
to  the  modest  selectivity  found  in  Ku  2:1  binding  is  only  observed  for  3-Man  at  25  °C  {K3  =  282 
M"'),  although  a  value  of  ^3  =  3  M"'  for  p-GIc  can  be  extrapolated  from  a  van’t  Hoff  plot. 
Interestingly,  a  second  receptor  association  does  not  occur  with  a-Man  suggesting  that  the  octyl 
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Fig.  6  a)  'H  NMR  (900  MHz,  CDCI3,  25°C)  spectra 
obtained  upon  the  titration  of  a  62.5  mM  solution  of  1 
to  a  0.98  mM  solution  of  p-Man  with  dashed  lines 
illustrating  the  induced  changes  in  5.  b)  Chemical 
shifts,  5,  of  p-Man  as  a  function  of  molar  equivalents 
of  1  at  5°C.  The  theoretical  global  fits  with  a  model 
incorporating  Kdimer,  K^,  and  K2  (dashed  line)  and  with 
Kdimer,  K^,  K2,  and  K3  (solid  line)  are  shown. 


chain  at  the  anomeric  position,  which  would  be  orientated  away  from  the  biphenyl  base  of  1  in 
1:a-Man  (Figure  5A),  interferes  with  the  association  of  a  second  receptor. 

The  cooperativity  that  facilitates  the  formation  of  l2:|3-Man  can  be  understood  through 
the  interaction  parameter,  a,  which  is  the  ratio  between  the  values  of  K3  in  the  presence  and  in 
the  absence  of  cooperativity,  the  latter  being  the  reference  Kg.  Since  the  two  faces  of  (3-Man  are 

inequivalent,  the  two  receptors  do  not  bind  to  identical 
sites,  thus  is  not  an  appropriate  reference  Kg.  Rather, 
since  K3  describes  the  binding  of  the  |3-face  of  (3-Man, 
which  contains  H\  H^,  and  H^,  in  the  presence  of  a 
receptor  bound  to  the  a-face,  which  contains  H'^,  the 
reference  Kg  would  describe  the  receptor  binding  to  the 
3-face  in  the  absence  of  a  receptor  bound  to  the  a-face. 
From  the  -63 'C  ^H-^FI  ROESY  spectrum  of  1:p-Man2, 
there  are  no  observable  NOEs  between  the  pyranoside 
protons  on  the  p-face  and  the  biphenyl  base  of  receptor 
1.  Likewise,  a  ROESY  spectrum  performed  at 
conditions  that  would  produce  significant  quantities  of 
the  1 :1  1:3-Man  complex  revealed  only  a  NOE  between 
Fl'^  of  the  a-face  and  H®.  Under  these  experimental 
conditions,  a  1:1  binding  event  between  1  and  the  p- 
face  does  not  occur  to  any  appreciable  extent,  so  an 
approximate  baseline  value  of  Kg  <  104  M'^  was 
obtained.  Thus,  an  estimated  K3  >  282  M"'  is  evidence 
for  positive  cooperativity  in  K3  with  p-Man,  with  a 
corresponding  a  of  at  least  2.7. 

To  determine  why  K3  occurs  exclusively  with  p- 
Man  at  room  temperature,  AH°  and  AS°  of  association 
for  each  binding  event  between  1  and  all  pyranosides 
were  obtained  by  repeating  the  NMR  titrations  at  20, 
15,  10,  and  5°C  and  subjecting  the  resulting  KgS  to  van’t  Hoff  analyses  (Table  3).  All 
monosaccharides  fit  satisfactorily  to  a  Kdimer  -1-  K^  model  at  all  temperatures,  except  for  p-Man  at 
all  temperatures  and  p-GIc,  at  15,  10,  and  5°C,  which  required  inclusion  of  the  2:1 
receptor:pyranoside  equilibrium  (K3)  to  achieve  satisfactory  fits  to  the  titration  data.  When 
comparing  the  thermodynamic  parameters  for  the  formation  of  each  1 :1  complex,  a  decrease  in 
AH°  occurs  with  a  decrease  in  magnitude  of  AS°  which  can  be  rationalized  within  the  context  of 
enthalpy-entropy  compensation.  Notably,  the  two  p-monosaccharides  with  the  highest  1:1 
binding  enthalpies,  p-Man  and  p-GIc,  are  the  only  pyranosides  that  participate  in  K3.  The  large 
difference  in  AH°  between  1:(3-Man  and  1:P-Glc,  -20.5  and  -16.5  kcal  mol"'  respectively, 
suggests  that  the  high  selectivity  in  K3  for  p-Man  is  the  result  of  1  forming  more  noncovalent 
contacts  with  p-Man  than  with  p-GIc,  resulting  in  a  more  preorganized  1:1  complex.  The  AH°s 
for  the  formation  of  l2:P-Glc  and  l2:P-Man,  -6.6  and  -11.0  kcal  mol  '  respectively,  suggest  that 
significantly  more  noncovalent  contacts  are  formed  in  the  latter,  indicating  that  the  selectivity  for 
p-Man  in  the  second  binding  event  is  enthalpically  driven. 

The  binding  geometries  of  1:p-Man  and  l2:P-Man  were  determined  and  compared  to 
understand  how  the  preorganization  of  1:P-Man  facilitates  the  formation  of  l2:P-Man.  A  'H-'H 
ROESY  experiment  that  was  performed  under  conditions  that  favor  the  1 :1  complex  (1  :p-Man  : 
1:p-Man2  =  2  :  1.5)  was  consistent  with  the  previously  determined  binding  geometry  for  1:p- 
Man.  Notably,  in  the  1:p-Man  geometry,  a  significant  portion  of  p-Man  is  exposed  to  solvent 
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Fig.  7  Energy  minimized  structure 
(AMBER*)  for  l2:p-Man  obtained  by  a 
mixed  low-mode/torsional  Monte  Carlo 
conformation  search.  The  carbons  of  1 
that  binds  the  a-face  are  colored  green, 
and  the  carbons  of  1  bound  to  the  p-face 
are  colored  red  for  clarity,  nitrogens  are 
blue,  oxygens  are  red,  and  intermolecular 
hydrogen  bonds  are  denoted  by  dashed 
lines.  Upon  the  association  of  the  second 
equivalent  of  1,  3  new  H-bonds  are 
observed  with  the  p-face. 


because  one  of  the  arms  of  1  is  orientated  downwards 
underneath  the  n-octyl  chain  of  P-Man,  which  provides 
a  window  for  a  second  equivalent  of  1  to  bind  onto  the 
exposed  p-face  of  P-Man.  The  structure  of  l2:P-Man 
was  determined  by  performing  a  ^H-^H  ROESY 
experiment  with  a  concentrated  2:1 
receptor:pyranoside  CDCI3  mixture  at  -10°C  where 
l2:P-Man  would  be  the  major  species  in  solution. 

NOEs  were  observed  between  protons  located  on  the 
biphenyl  base  and  pyrroles,  suggesting  a  geometry 
where  one  molecule  of  p-Man  is  encaged  by  two 
molecules  of  1,  and  the  two  molecules  of  1  are  in 
close  contact.  By  using  a  similar  AMBER* 
conformational  search  that  was  used  to  obtain  the  1  :P- 
Man  structure,  a  l2:P-Man  structure  was  obtained 
(Figure  7)  that  is  consistent  with  the  thermodynamic 
data  and  the  ID  and  2D  NMR  spectra.  In  contrast, 
the  ^H-^H  ROESY  of  p-GIc  and  1  under  conditions  that 
favor  the  1 :1  complex  revealed  that  the  axial  protons 
on  both  faces  of  p-GIc  have  observable  NOEs  with  the 
biphenyl  base  of  1,  indicating  that  more  than  one 
binding  geometry  exists  at  equilibrium  since  1  cannot 
interact  with  both  faces  of  p-GIc  simultaneously. 

Because  the  entropic  penalty  for  the  preorganization 
of  1:p-Glc  has  not  been  paid,  we  speculate  that  the 
AS°  associated  with  a  second  molecule  of  1  binding  to 
the  1:1  complex  becomes  prohibitively  high, 
preventing  the  formation  of  a  2:1  complex.  It  should  be 
noted  that  in  the  calculated  structure  of  l2:P-Man  only 
three  arms  of  each  equivalent  of  1  participate  in 
binding  to  p-Man,  so  the  fourth  arms  can  be  replaced 
in  future  receptors  to  increase  function. 

Implications  of  competing  equilibria  on  the  selectivity  of  1 

In  the  presence  of  p-Man,  the  1:1  complex,  1:p-Man,  can  either  bind  a  second  molecule  of 
pyranoside,  to  form  1:p-Man2  (K'2),  or  1,  to  form  l2:P-Man  (K'3).  The  outcome  of  these  two 
competing  processes  is  dependent  on  the  concentrations  of  each  substrate,  with  K2  dominating 
at  high  concentrations  of  p-Man  relative  to  1,  whereas  K'3  is  preferred  the  presence  of  excess  1. 
Thus,  the  numerous  equilibria  operating  simultaneously  in  a  solution  of  1  and  pyranoside  (K'dimer, 
Ku  K'2  and  K'3)  and  the  concentration  dependence  of  cooperative  binding  complicate  the  analysis 
of  the  overall  selectivity  of  1  for  the  eight  pyranosides.  The  BC50  values  of  each  of  the 
pyranosides  -  defined  as  the  total  concentration  of  receptor  1  needed  to  bind  50%  of  the 
available  pyranosides  in  solution  -  were  computed  over  a  pyranoside  concentration  range  of  0 
to  10.0  mM  to  probe  both  the  magnitude  and  concentration  dependence  of  the  selectivity  of  1 
for  the  eight  pyranosides  (Figure  8A).  Note  that  all  equilibria  present  in  a  receptor/pyranoside 
mixture  contribute  toward  the  calculated  BC50  value,  and  a  higher  I/BC50  value  is  indicative  of  a 
higher  overall  binding  strength. 


Fig.  8  I/6C50  plots  describing  a)  the  pyranoside 
concentration  dependence  for  the  overall 
binding  strength  (I/6C50)  of  receptor  1  toward 
each  pyranoside,  where  higher  values  are 
indicative  of  stronger  binding  and  crossover 
points,  marked  within  a  circle  for  p-GIc^a-Man 
and  square  for  p-GIc^p-Man,  denote  changes 
in  binding  preference  with  increasing 
pyranoside  concentration;  and  b)  the  influence 
of  each  positive  cooperative  binding  equilibria 
(K2  and  K'3)  on  the  BCso  value  of  receptor  1  for 
3-Man. 
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The  I/BC50  plots  reveal  the  effects  of  positive  cooperativity  on  the  binding  affinity  of  1 
toward  mannosides.  Molecules  with  K2  possess  parabolic  curves  in  the  I/BC50  plot.  This  is 
particularly  evident  for  a-Man  which  possesses  the  highest  affinity  for  1  compared  to  all  other 
pyranosides  at  a  0.8  mM  pyranoside  concentration.  Interestingly,  for  (3-Man,  which  exhibits 
cooperativity  in  both  K2  and  K3,  a  broader  parabolic  shape  is  observed,  and  little  concentration 
dependence  on  the  overall  binding  strength  occurs  between  0  and  4.0  mM  of  3-Man.  When  the 
calculated  I/BC50  plot  is  compared  to  a  hypothetical  plot  in  which  only  K^  and  K3  are  present,  it 
was  found  that  the  presence  of  K3  induces  a  substantial  increase  in  selectivity  at  low 
concentrations  (<  1.5  mM)  of  3-Man  (Figure  8B).  Alternatively,  when  a  hypothetical  plot  is 
generated  with  only  and  K2,  K2  enhanced  selectivity  when  [(3-Man]  is  greater  than  1.0  mM, 
thus  demonstrating  the  direct  contribution  of  cooperative,  compiex  equiiibria  on  seiectivity  in 
saccharide  receptors  (Figure  8B).  Moreover,  these  plots  confirm  that  pyranoside  concentration 
controls  which  equilibrium,  i.e.  K2  or  K3,  prevails.  At  concentrations  below  1.0  mM  of  3-Man, 
higher  stoichiometry  binding  is  achieved  through  K3,  while  at  greater  3-Man  concentrations,  K2 
dominates. 

As  each  I/BC50  plot  is  a  measure  of  affinity  between  1  and  pyranosides,  comparing  the 
plots  of  two  or  more  pyranosides  provides  a  means  of  assessing  selectivity.  Thus,  the 
concentration  of  a-Man  where  the  maximum  I/BC50  value  is  obtained  (0.8  mM)  is  also  where 
the  maximum  selectivity  occurs  for  a-Man  relative  to  the  other  pyranosides.  Beyond  this  value, 
the  selectivities  gradually  approach  2:1  mannoside:pyranoside,  which  is  a  consequence  of  the 
differences  in  binding  stoichiometry  -  receptor  1  can  accommodate  two  molecules  of 
mannosides  and  only  binds  1 :1  for  the  other  pyranosides.  Importantly,  crossover  points  in  the 
BC50  plots  are  observed  between  3-Glc  and  the  mannosides  (Figure  8A),  meaning  that  the 
seiectivity  of  1  changes  as  a  function  of  pyranoside  concentration,  which  is  a  form  of  moiecuiar 
iogic,  and  a  major  goai  of  Y3.  At  low  pyranoside  concentrations  (<  0.3  mM),  receptor  1  binds  3- 
Glc  with  the  highest  affinity  owing  to  its  high  value  of  Ku  which  effectively  competes  with  the 
cooperative  binding  of  1  toward  mannosides.  From  0.3  mM  to  1 .0  mM,  the  binding  order  for  1  is 
a-Man  >  3-Glc  >  3-Man.  At  higher  pyranoside  concentrations  (>  1 .0  mM),  the  binding  order  for  1 
changes  to  a-Man  >  3-Man  >  3-Glc  as  a  result  of  the  2:1  binding  stoichiometry  that  occurs  only 
with  mannosides.  Consequently,  a  new  concentration  dependence  has  been  discovered  that 
governs  the  stoichiometry  of  the  resulting  complex  and  the  pyranoside  preference  of  the 
receptor. 

Significance  and  outiook 

The  flexible  receptor  1  achieves  selectivity  between  pyranosides  that  may  differ  only  by  the 
orientation  of  a  single  hydroxyl  group,  despite  the  entropic  penalty  that  must  be  paid  to  organize 
the  complexes.  In  fact,  the  lack  of  preorganization  in  the  strong  1:3-Glc  complex  precludes  the 
formation  of  a  stable  12:3-Glc  structure,  so  entropy  contributes  significantly  towards  selectivity. 
Although  receptor  1  does  not  achieve  the  same  overall  affinity  as  the  best  rigidly  preorganized 
mannose-specific  receptors,  the  selectivity  is  comparable  despite  the  differences  in  binding 
mechanisms,  thus  confirming  the  hypothesis  that  increasing  the  receptor  dynamics  reveals  new 
binding  geometries  because  of  the  ability  of  flexible  hosts  to  dynamically  explore  conformational 
space.  Synthetic  carbohydrate  receptors  can  provide  insight  into  the  subtleties  of  natural  lectin- 
carbohydrate  interactions  and  reveal  how  carbohydrate  recognition  conveys  complex 
information  in  biological  networks.  Eukaryotic  cell  surfaces  are  coated  with  a  carbohydrate 
layer,  the  glycocalyx,  where  the  multivalent  presentation  of  carbohydrates  on  a  cell  surface 
enhances  binding  affinity,  a  phenomenon  termed  the  cluster  glycoside  effect.  Interestingly,  the 
selectivity  of  1  changes  with  pyranoside  concentration,  where  1  preferentially  binds  glucosides 
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at  low  concentration  (<0.3  mM)  and  mannosides  at  higher  concentrations  that  are  similar  to  the 
increased  localized  concentrations  present  in  mannoside  clusters  on  cell  surfaces.  While 
concentration  dependent  selectivity  may  be  unprecedented  with  synthetic  carbohydrate 
receptors,  multivalency  and  cooperativity  are  ubiquitous  in  biology,  so  concentration  dependent 

switching  may  commonly  occur  with  natural 
lectins  and  could  have  implications  for 
hierarchical  organization  and  information 
transfer  in  biological  networks.  This  work 
{Chemical  Science,  2013,  4,  357-367) 

represents  the  accomplishment  of  the  major 
Y1  and  Y2  goals  with  respect  to  synthetic 
lectins  by  creating  a  new  receptor  that 
engages  cooperativity,  a  binding  mode 
ubiquitous  to  natural  lectins,  to  perform 
molecular  logic  operations. 

C.2.  Towards  Graphene  Based  Glycan 
Arrays 

In  eukaryotic  cells,  the  primary  location  of 
carbohydrate  recognition  is  on  the  cell  surface 
glycan  layer,  termed  the  glycocalyx.  To 
understand  how  surface  properties,  such  as 
density  and  orientation,  affect  carbohydrate 
recognition,  we  are  exploring  new  chemistries 
and  surfaces  to  immobilize  carbohydrates  into 
patterns  with  controlled  density  and 
orientation.  The  first  of  these  surfaces  is 
single-layer  graphene  (SLG),  which  is 
composed  of  sheets  of  all  sp^-hybridized 
carbon.  The  recent  report  by  Haddon  et  al. 
demonstrating  that  SLG  participates  in  Diels- 
Alder  reactions  as  a  dienophile  at 
temperatures  as  low  as  50  °C  over  3  h 
suggested  to  us  conditions  that  could  be  used  to  scalably  pattern  graphene  at  ambient 
temperatures  and  atmosphere.  Because  of  their  negative  activation  volumes,  cycloaddition 
reactions  are  significantly  accelerated  in  pressurized  reaction  vessels,  thus,  we  reasoned  that  a 
localized  force  between  SLG  and  a  diene  (Figure  9a)  would  accelerate  the  Diels-Alder  reactions 
and  thereby  immobilize  molecules  site-specifically  onto  the  basal  plane  of  graphene  through  the 
formation  of  two  new  C-C  bonds  (Figure  9b). 


SLG  2 


Fig.  9  Tip-induced  force  accelerated  Diels-Alder 
reaction  on  a  single  layer  graphene  (SLG)  surface,  a) 
An  elastomeric  tip-array,  b)  The  tip-array  coated  with 
an  ink  mixture  (red)  consisting  of  a  cyclopentadiene 
and  poly(ethylene  glycol)  (PEG),  c)  The  inked  tip  array 
is  pushed  into  the  SLG  surface,  d)  Following  rinsing  of 
the  surface  to  remove  the  PEG  and  excess 
cyclopentadiene,  only  the  covalently  immobilized 
molecules  remain  on  the  surface,  e)  Diels-Alder 
reaction  between  functionalized  cyclopentadiene  and 
an  SLG  surface,  f)  Cy3-containing  Raman  active  1  and 
ferrocene-containing  electrochemically  active  2  ink 
molecules  used  to  confirm  force-accelerated 
patterning. 
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To  demonstrate  that  force  accelerated  cycloadditions  could  covalently  pattern  large 
areas  (~1  cm^)  of  SLG  sheets,  we  used  an  elastomeric  tip  array  mounted  onto  the  piezoelectric 
actuators  of  an  atomic  force  microscope  (AFM)  to  site- 
specifically  apply  a  force  between  functionalized 
cyclopentadienes  and  SLG  sheets.  These  tip  arrays  are 
commonly  used  for  Polymer  Pen  Lithography,  where 
patterns  are  formed  by  ink  transfer  from  the  tips  to  the 
surface  through  an  aqueous  meniscus,  resulting  in  a 
relationship  between  dwell  time  and  feature  size. 

Moreover,  their  large  areas  (>1  cm^)  and  the  computer 
controlled  movement  of  the  piezoactuators  that  hold  the 
array  provide  high  throughput  and  flexible  pattern 
designecause  the  elastomeric  tips  also  compress  upon 
contact  with  surfaces,  they  can  apply  a  predictable 
force  between  molecular  inks  and  a  surface.  The 
relationship  between  applied  force  and  the  resulting 
feature  sizes  (Eq.  1), 


F  = 

where  L 


{Lf-Lt)E2Lt 


(Eq. 1) 

1/ 

.f  is  feature  edge  length  of  the  top  of  the  tip, 
is  the  feature  edge  length  of  the  bottom  of  the  tip,  E2  is 
the  compression  modulus  of  poly(dimethyl  siloxane) 

(PDMS),  and  v  is  the  Poisson’s  ratio  of  PDMS,  can  be 
used  to  determine  the  force  between  the  tips  and  the 
surfaces.  As  a  result,  in  this  experimental  design,  the 
position,  force,  and  time  can  all  be  controlled  precisely 
to  pattern  surfaces  with  micrometer  scale  features  over 
a  square  centimeter  area. 

Raman-active  cyanine  3  (Cy3)  containing 

cyclopentadiene  1  and  electrochemically-active 
ferrocene  cyclopentadienes  2  (Fig.  9c)  were  designed 
to  characterize  the  bonding  and  density  upon  reaction 
between  the  SLG  surface  and  the  cyclopentadienes. 

Cyclopentadienes  react  quickly  in  Diels-Alder  reactions  compared  to  open  chain  dienes 
because  they  are  structurally  preorganized  for  reaction,  and  as  a  result,  they  have  been  utilized 
already  in  the  context  of  surface  patterning.  Raman  spectroscopy  has  become  the  standard  tool 
for  characterizing  chemical  modifications  onto  the  basal  plane  of  graphene.  The  Fladdon  group 
observed  that  following  a  Diels-Alder  reaction  onto  graphene,  the  D  band  at  1345  cm'^  that 
corresponds  to  the  Aig  breathing  vibration  of  sp^  carbon  rings,  which  is  suppressed  in  pure 
graphene,  increases  significantly  because  the  introduction  of  defects  or  covalently  adsorbed 
molecules  reduces  the  symmetry  of  the  graphene  lattice.  As  a  result,  the  ratio  of  the  D-  and  G- 
band  intensities  (/d//g)>  determined  using  the  integration  of  the  peaks,  is  a  relative  measure  of 
the  degree  of  functionalization  of  graphene.  Alternatively,  electrochemistry  can  confirm  the 
immobilization  of  the  ink  onto  the  surface  and  quantify  the  density  of  surface-bound  molecules. 


Fig.  10  a)  Light  microscopy  image  (10x 
magnification)  of  2  x  3  dot  arrays  of  a 
mixture  containing  1  and  PEG  with  varying 
dwell  times  (30,  27,  24,  21,  18,  15  min) 
were  patterned  by  each  pen  in  the  PPL  tip 
array.  Scale  bar  is  200  pm.  Inset  is  a 
magnified  image  of  one  array.  Scale  bar  is 
20  pm.  b)  Raman  mapping  image  (1324 
cm”\  D  band)  of  2  x  3  dot  arrays  of  1 
covalently  immobilized  onto  the  SLG. 
Scale  bar  is  20  pm. 
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Molecules  1  and  2  were  synthesized  and 
characterized  by  NMR,  NMR,  and  high- 
resolution  mass  spectrometry,  and  all  analytical  data 
were  consistent  with  the  proposed  structures.  The 
8500  tip  arrays  with  a  tip-to-tip  spacing  of  80  or  160  pm 
were  prepared  following  previously  published  literature 
protocols  and  are  composed  of  PDMS  pyramids 
mounted  onto  a  glass  support.  To  induce  the  Diels- 
Alder  reaction  between  1  and  the  SLG  surface,  1  (0.8 
mg,  1.2  mmol)  and  poly(ethylene  glycol)  (PEG)  (2000  g 
mor\  10  mg  ml'^)  in  0.8  mL  60  :  20  THF:H20,  which 
was  sonicated  to  ensure  solution  homogeneity,  were 
spin  coated  (2000  rpm,  2  min)  onto  a  tip  array.  The 
PEG  matrix  that  encapsulates  the  cyclopentadienes 
ensures  even  distribution  across  the  tip  array,  and  in 
the  case  of  Polymer  Pen  Lithography,  transport  from 
the  tips  to  a  surface  is  predictable  and  reproducible. 
The  tips  were  then  mounted  onto  the  z-piezo  of  an 
AFM  that  was  specially  equipped  with  an  apparatus  to 
hold  the  tip  arrays,  an  environmental  chamber  to 
regulate  the  humidity,  and  customized  lithography 
software  to  control  the  position,  force,  and  dwell-time  of 
the  tips.  A  2x3  dot  pattern  of  1  with  feature-to-feature 
spacing  of  20  pm  was  produced  by  each  tip  in  the  array 
by  pushing  the  tips  into  the  SLG  surface  (SLG  on  285 
nm  Si02)  at  times  ranging  from  15-30  min  and  a  force 
of  -100  mN  at  each  spot.  The  transfer  of  the  1/PEG 
mixture  to  the  surface  was  confirmed  by  light 
microscopy  (Figure  10a),  and  the  2x3  1/PEG  patterns 
and  the  approach  dot  used  to  level  the  tip  array  with 
respect  to  the  surface  are  clearly  visible. 

After  washing  the  surfaces  with  EtOH  and  H2O 
to  remove  unreacted  1  and  PEG,  the  surface  bonding 
was  analyzed  by  Raman  microscopy  (Renishaw  inVia, 
633  nm  excitation).  A  Raman  map  of  the  surface  that 
was  obtained  following  force  accelerated  printing  of  1 
revealed  a  2x3  pattern  of  features  where  Id  was  elevated  significantly  compared  to  surrounding 
areas  (Figure  10b).  The  dimensions  and  feature  sizes  of  these  2x3  pattern,  with  20  pm  spacing 
between  features,  matched  with  the  pattern  of  features  printed  by  the  pen  array.  The  elevated  Id 
was  observed  at  all  points  where  the  tips  were  pressed  into  the  surface  for  all  dwell  times. 
Importantly,  control  experiments  where  1  was  not  present  in  the  ink  mixture  or  where  1  was 
present,  but  force  was  not  applied  to  the  surface  upon  ink  transfer,  did  not  produce  similar 
patterns  or  significantly  elevated  IdUg  in  the  Raman  maps,  confirming  that  the  diene  is 
necessary  for  changes  in  bonding  to  occur.  The  Raman  spectra  associated  with  different  points 
on  this  map  further  confirmed  that  the  changes  in  bonding  were  produced  because  of  the 
occurrence  of  localized  Diels-Alder  reactions  (Figure  11).  A  Raman  spectrum  taken  at  a  point 
where  the  tips  were  pressed  into  the  surface  had  peaks  corresponding  to  the  SLG  D-band  and 
1,  as  well  as  an  increased  /d//g  value  of  0.95,  compared  to  0.49  for  the  unaltered  surface,  0.46 
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Fig.  11  a)  Raman  spectrum  of  1  on  an 
Si02  surface,  b)  Raman  spectrum  from  a 
map  with  force  accelerated  printing  of  1 
(Fig.  2b)  taken  at  a  position  with  increased 
D-band  intensity  (/□)  compared  to 
unreacted  SLG  also  has  peaks 
corresponding  to  1.  c)  Raman  spectrum 
from  the  map  of  the  same  printed  surface 
taken  at  a  position  without  increased  fc.  d) 
Raman  spectrum  of  control  experiment 
where  PEG  and  force  were  applied  to  the 
SLG  surface  without  1,  taken  at  a  point 
where  the  tips  were  pushed  into  the  SLG 
surface,  e)  Raman  spectrum  of  pure  SLG. 
Red  lines  mark  peaks  of  1  (11124,  1157, 
1178,  1272,  1387,  1594  cm’’').  Blue  lines 
mark  SLG  peaks  (1345,  1584  cm“Y 


where  the  tips  had  been  pressed  into  the  surface  in  the  absence  of  1 ,  and  0.39  for  the  original 
SLG  surface.  The  changes  in  the  IqIIg  in  the  Raman  maps  indicate  1  is  on  the  surface  and 
changes  in  bonding  from  sp^  to  sp^  that  are  consistent  with  those  previously  observed  for  Diels- 
Alder  reactions  on  SLG  surfaces,  occur  in  the  same  position.  Additionally,  control  experiments 
confirm  that  these  changes  in  the  spectrum  only  happen  under  conditions  where  the  Diels-Alder 
reaction  can  proceed.  Unlike  cycloadditions  under  pressure,  where  rate  accelerations  arise 
because  of  the  negative  activation  volume,  it  is  conceivable  that  rate  acceleration  may  also 
arise  because  of  the  distortions  of  the  SLG  rr-bonds  under  applied  force  that  increase  their 
reactivity. 

Electrochemically  active  cyclopentadiene  2  was  patterned  onto  SLG  following  a  similar 
protocol  described  above,  and  the  immobilization  density  of  2  on  the  SLG  surface  was  analyzed 
by  cyclic  voltammetry  (CV).  Each  tip  in  the  pen  array  produced  a  2  x  3  dot  pattern  over  the  1 
cm^  area  covered  by  the  tip  array  with  a  dwell  time  of  30  s  at  each  spot,  and  ink  deposition  was 

confirmed  by  optical  microscopy  with  an  average 
feature  edge  length  of  7.1  pm  and  area  of  50.4 
pm^.  Following  washing  of  the  surface  with  EtOH 
and  H2O  to  remove  the  PEG  and  unreacted  2, 
CV  was  carried  using  a  Ag/AgCI  reference 
electrode,  a  Pt  counter  electrode,  and  the 
patterned  SLG  as  the  working  electrode.  A  strong 
redox  peak  at  £°  =  590  mV  (vs.  Ag/AgCI) 
confirmed  the  presence  of  the  ferrocene 
(fc)/ferrocenium  (fc"")  reversible  redox  couple 
from  2  (Figure  12),  which  is  shifted  anodically 
compared  to  fc  because  of  the  ester  linking  the  fc 
to  the  cyclopentadiene.  The  linear  relationship 
between  peak  current  and  scan  rate  confirmed 
that  2  is  immobilized  on  the  SLG  surface,  but  that 
the  localized  changes  in  bonding  from  sp^  to  sp^ 
do  not  prevent  conduction  through  the  SLG. 
However,  the  difference  between  oxidation  and 
reduction  peaks,  AE,  may  indicate  an  increase  in 
resistance  upon  changes  in  chemical  bonding  from  sp^  to  sp^.  The  surface  density  of  fc  within 
each  feature,  T ^0,  was  determined  from  the  CV  measurements.  A  Tfc  of  (5.34  ±  0.76)  x  10^'*  cm'^ 
was  obtained  and  if  the  density  of  Ti-bonds  on  the  graphene  surface  is  considered,  this  value  of 
r ^0  corresponds  to  approximately  29%  bond  functionalization.  Control  experiments  where  2  was 
deposited  without  force  did  not  result  in  any  observable  current  corresponding  to  the  fc/fc"^  redox 
couple  after  washing,  confirming  that  force  is  necessary  to  induce  the  Diels  Alder  reaction  under 
these  conditions.  In  conclusion,  SLG  sheets  on  Si02  substrates  have  been  patterned  covalently 
with  organic  small  molecules  through  a  force-accelerated  Diels-Alder  reaction  induced  by  an 
array  of  pyramidal  elastomeric  tips.  The  results  of  these  efforts  are  described  in  further  detail  in 
an  article  recently  published  in  the  Journal  of  the  American  Chemical  Society  (J.  Am.  Chem. 
Soc.  2013,  135,  9240). 
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Fig.  12  Cyclic  voltammetry  (CV)  characterization 
of  the  SLG  patterned  with  2  using  a  Pt  counter 
electrode  and  Ag/AgCI  reference  electrode  in  0.1 
M  HCIO4  electrolyte  solution.  The  inset  shows  the 
linear  relationship  between  scan  rate  and  current, 
indicating  that  2  is  surface  immobilized. 
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C.3.  Spatially  controlled 

photochemistry,  and  its  utilization  in 
the  synthesis  of  multivalent  glycan 
arrays 

Ideally  one  would  like  to  generate 
features  that  behave  as  multivalent 
scaffolds,  that  can  be  prepared  from 
readily  available  starting  materials,  and 
that  can  be  patterned  using  easy  to 
implement  surface  chemistry.  Herein,  we 
describe  a  new  approach  that  combines 
PPL  or  BPL  with  acrylate  or  methacrylate 
photopolymerization  chemistry  to 
generate  arrays  of  microscale  features  of 
brush  polymers  that  are  side-chain 
functionalized  with  fluorophores  or 
carbohydrates.  These  polymers  exhibit 
excellent  GBP  binding  ability  by 
mimicking  the  multivalent  glycan 
structures  typically  found  on  cell 
surfaces.  Using  the  thousands  of  tips  in  a 
BPL  array  to  localize  light  on  a  surface, 
the  exposure  time  at  each  position  on  the 
surface  during  the  photopolymerization 
could  be  varied  precisely,  providing  in 
situ  control  over  polymer  length  during 
patterning.  Thus,  in  addition  to  achieving 
orders  of  magnitude  improvement  over 
binding  sensitivity  as  a  result  of 
multivalency,  the  glycan  arrays  prepared 
by  this  3D  photolithography  demonstrate 
that  a  relationship  exists  between  polymer  length  and  GBP  limit  of  detection. 

Inspired  by  the  thiol-ene  reaction,  which  has  been  used  successfully  in  the  context  of 
carbohydrate  immobilization,  we  turned  to  the  thiol-acrylate  photopolymerization  as  a  potential 
route  towards  creating  glycan-functionalized  polymer  brushes  with  pendant  carbohydrate  side 
groups  because  of  (1)  the  compatibility  of  the  radical  reaction  mechanism  with  the  functional 
groups  common  to  carbohydrates,  and  (2)  the  ability  to  induce  the  reaction  with  light  made  it  an 
appropriate  choice  for  high-throughput  processing  with  BPL.  In  the  presence  of  the  photoinitiator 
2,2-dimethoxy-2-phenylacetophenone  (DMPA),  UV  light  induces  photopolymerization  of  acrylate 
or  methacrylate  monomers  from  a  thiol  initiator.  By  combining  this  polymerization  with  BPL 
arrays,  control  over  feature  diameter,  feature  height  and  feature  spacing  was  possible. 
Moreover,  the  ease  with  which  acrylate  and  methacrylatemonomers  are  functionalized  suggests 
that  this  approach  could  find  broad  use  beyond  glycan  arrays.  However,  the  photoinduced  thiol- 
initiated  acrylate  polymerization  has  not  been  explored  extensively  as  a  tool  to  synthesize  brush 
polymers,  and  there  are  no  examples  of  glycan  brush  polymers  prepared  by  the  thiol-initiated 
acrylate  radical  photopolymerization.  Therefore  before  glycan  arrays  could  be  synthesized,  it 
was  necessary  to  first  develop  methods  to  pattern  brush  polymers  by  combining  the  acrylate 
photopolymerization  reaction  with  BPL. 


Fig.  13  (a)  Process  for  inducing  the  thiol-ene  and  (meth)acrylate 
polymerization  reactions  on  a  surface  by  beam  pen  lithography 
(BPL).  i)  The  tip-array,  coated  with  an  ink  mixture  (blue)  containing 
the  probe  molecules  and  DMPA  in  a  poly(ethylene  glycol)  (PEG) 
matrix,  deposits  the  ink  mixture  onto  thiol-terminated  glass  slides,  ii) 
Light  passing  through  the  tips  of  the  beam  pens  locally  exposes  the 
patterned  surface  to  induce  hemolytic  cleavage  of  the  DMPA  and 
initiate  the  polymerization  of  the  monomers  within  the  PEG  matrix, 
iii)  Following  rinsing  of  the  surface  to  remove  the  PEG  and  excess 
molecules,  only  the  covalently  immobilized  polymers  remain. 
Methacrylate  brush  polymers  are  shown,  (b)  Molecules  prepared 
for  the  BPL  and  polymer  pen  lithography  (PPL)-induced  thiol-ene 
and  thiol-(meth)acrylate  polymerization  reactions. 
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The  BPL  induced  radical  photopolymerization  is 
a  new  type  of  3D  scanning  probe  nanolithography, 
where  the  feature  diameter,  feature  height,  feature 
position  can  all  be  controlled  independently  during 
printing.  The  photopolymerization  is  initiated  by 
photochemical  cleavage  of  the  DMPA  to  form  the  thiol 
radical,  which  in  turn  initiates  acrylate  propagation. 
These  thiol  radicals,  however,  lack  the  stability  of 
conventional  radical  initiators,  and  because  the  chain 
grows  from  the  thiol,  polymerization  stops  when  the 
light  is  turned  off.  Because  polymerization  proceeds 
only  during  illumination,  we  reasoned  that  we  could  use 
this  property  of  the  thiol-acrylate  photopolymerization 
to  control  brush  polymer  length  by  varying  the  dose  of 
light  delivered  to  each  patterned  feature.  In  the  current 
approach,  feature  position  is  controlled  by  the 
piezoactuator  that  moves  the  tip  arrays,  diameter  is 
varied  by  the  dwell-time  between  the  pen  arrays  and 
the  surface,  and  feature  height  is  a  function  of  the 
irradiation  time.  To  demonstrate  the  control  over  feature 
diameter  during  thiol-ene  and  radical  polymerization 
printing,  fluorescent  molecules  1  and  2  were  patterned 
into  a  gradient  by  PPL  with  dwell  times  of  50-  100  000 
ms  and  subsequently  irradiated  for  3  h  (Figure  13). 
Chain  growth,  in  combination  with  the  inability  of  thiol- 
initiators  to  sustain  polymerization  in  the  absence  of 
continuous  irradiation,  is  essential  for  obtaining  precise 
control  over  polymer  chain  length  with  exposure  time. 
Although  chain  transfer  is  likely  to  occur,  chain  transfer 
does  not  affect  the  growth  kinetics  of  brush  polymers. 
Thus,  3D  lithography  capable  of  immobilizing  polymers 
laden  with  functional  side-chains  -  while  directing 
feature  diameter,  feature  height,  and  feature  position 
with  sub-micrometer  resolution  -  is  achieved  by 
combining  BPL  with  the  photoinduced  radical 
polymerization.  Importantly,  this  is  the  first 
demonstration  of  patterning  by  depositing  material  with 
a  BPL  array  followed  by  exposure,  and  the  first 
example  of  covalent  photochemistry  with  BPL,  showing 
that  coordinated  material  deposition  and 
photochemistry  is  possible  through  this  technique. 

Our  primary  aims  in  developing  this  new  scanning  probe  lithography  were  to  create  a 
method  to  prepare  glycan  arrays  with  multivalent  carbohydrates  and  to  investigate  how 
molecular  architecture  affects  the  binding  affinity  between  patterned  glycans  and  GBPs.  We 
expected,  as  an  important  validation  of  this  new  approach,  that  the  high  glycan  valency  of  our 
brush  polymers  would  increase  significantly  the  binding  towards  ConA  compared  to  arrays 
composed  of  glycan  monolayers.  One  additional  advantage  of  this  new  approach  is  that  the 
pattern  versatility  enabled  by  combining  BPL  with  in  situ  glycopolymer  synthesis  can  be  used  to 
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Fig.  14.  (a)  Fluorescence  image  (Nikon  Eclipse 
Ti,  Aex=  532-587  nm,  Aem=  608-687  nm)  of  an 
array  of  1  printed  by  BPL  with  illumination  times 
of  2,  5,  10,  and  20  min.  Inset  shows  an  array 
prepared  by  a  single  tip.  (b)  Intensity  profile  of 
the  features  printed  with  different  illumination 
times  indicated  by  the  white  line  in  (a).  Inset 
shows  the  cartoon  image  of  the  polymers 
formed  on  the  surface,  (c)  The  relationship 
between  exposure  time  and  fluorescence 
intensity  or  height.  The  variation  in  measured 
height  was  5%,  1%,  7%  and  43%  from  2  in  to 
20  min,  respectively. 


produce  multiplexed  arrays  of  polymers  of 
systematically  varied  lengths  on  the  same 
surface.  Glucosides  3  and  4  and  mannoside  5 
were  printed  by  either  the  thiol-ene  or  acrylate 
polymerization  reactions  onto  thiol-terminated 
glass  slides,  and  the  resulting  arrays  were 
assayed  against  solutions  containing  different 
concentrations  of  the  GBP  concanavilin  A 
(ConA).  ConA  was  selected  for  these  studies 
because  it  possesses  many  of  the  characteristic 
binding  properties  of  GBPs  that  complicate  the 
analysis  of  GBP-glycan  interactions,  including 
(1)  multiple  recognition  sites,  (2)  cross 
specificity,  and  (3)  low  1  :  1  binding  affinity. 

Specifically  ConA  possesses  two  carbohydrate 
recognition  domains  separated  by  a  distance  of 
6.8  nm  and  binds  both  a-mannosides  (Ka  10^- 
10"^  M- )  and  a-glucosides  (Ka  10^  M-),  albeit 
with  a  preference  for  a-mannosides  (Figure  14). 

The  signal  for  the  patterns  of  the 
multivalent  glucoside  3  that  were  prepared  by 
the  acrylate  photopolymerization  are 
approximately  20  times  greater  than  those 
obtained  with  amannoside  monolayers  by  the 
thiol-ene  and  the  CuAAC  reactions.  Arrays  of  3 
prepared  by  the  BPL-induced  thiol-acrylate 
photopolymerization  decreased  steadily  with 
decreasing  Cy3-ConA  concentration,  but  still 
had  appreciable  signal  at  0.43  mM,  which  is  at 
least  one  order  of  magnitude  more  sensitive 
than  the  glycan  monolayers  prepared  by  the 
thiol-ene  reaction,  and  nearly  two  orders  of 
magnitude  more  sensitive  than  a-mannoside 
monolayers  prepared  by  the  CuAAC  click 
reaction  (Figure  15).  The  role  of  chain  length  on 
binding  was  explored  by  determining  the 
detection  threshold  of  the  different  brush  polymers  of  3,  whose  lengths  were  varied  by  altering 
the  UV  irradiation  time.  The  feature  heights  of  the  ConA-brush  polymer  aggregates  increase 
approximately  linearly  with  exposure  time,  further  supporting  the  chain  growth  polymerization 
mechanism.  This  work  validates  three  important  concepts:  (1)  the  same  tips  can  be  used  to 
pattern  materials  and  direct  light  to  them  to  drive  photochemistry,  (2)  combinatorial  covalent 
photochemistry  can  be  performed  using  BPL,  and  (3)  arrays  of  biomaterials  constructed  through 
these  approaches  can  be  used  to  address  challenging  scientific  problems,  in  this  case 
understanding  GBP  binding.  The  results  of  these  efforts  are  described  in  further  detail  in  an 
article  recently  published  in  Chemical  Science  [Chem.  Sci.  2014,  5,  2023). 


Fig.  15  (a)  Fluorescence  image  (Nikon  Eclipse  Ti,  Aex= 
532-587  nm,  Aem=  608-687  nm)  of  an  array  of  3  printed  by 
BPL  with  illumination  times  of  2,  5,  10,  and  20  min. 
Following  washing,  the  surfaces  were  exposed  to  a  21 .7 
pM  solution  of  Cy3-labeled  ConA.  Inset  shows  an  array 
prepared  by  a  single  tip.  (b)  Binding  of  Cy3-ConA  to 
glycan  arrays  prepared  using  different  surface  chemistries 
and  technologies,  including  a-mannose  immobilized  by 
the  CuAAC  reaction  by  PPL  (green),  a-glucose  by  the 
acrylate  polymerization  of  3  by  BPL  (red),  thiol-ene  of  a- 
mannose  alkene  by  PPL  (blue),  thol-ene  reaction  with  4 
by  PPL  (purple),  and  binding  of  the  acrylate  surface  of  3 
with  PNA  (black)  by  PPL. 
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C.4.  Role  of  Force  on  the  Kinetics  of  the  Hi 

The  mechanisms  by  which  force 
accelerates  surface  reactions  remains  poorly 
understood,  because  there  are  few 
experimental  methods  capable  of  controlling 
precisely  critical  reaction  parameters  such  as 
reaction  time  t  and  applied  P.  Forces  can  be 
controlled  accurately  using  a  single  AFM  tip, 
but  the  low  throughput  of  this  method 
necessitates  complex  single-molecule 
analyses.  Alternatively  soft  lithography  has 
been  used  to  investigate  force-accelerated 
reactions  over  large  areas,  which  simplifies 
analysis,  but  the  applied  P  is  difficult  to 
control.  As  a  result,  these  AFM  and  soft 
lithography  experiments  can  lead  to 
inconsistent  and  contradictory  results.  For 
example,  attempts  to  use  force  to  induce  the 
Huisgen  cycloaddition  reaction  between 
alkyne-terminated  surfaces  and  azide 
substrates,  or  vice  versa,  in  the  absence  of 
Cu'  have  alternately  concluded  that  the 
reaction  does  and  does  not  proceed  under 
applied  pressure.  Before  the  potential  of 
force-induced  reactions  on  surfaces  can  be  realized,  there  is  a  need  for  new  experimental 
methods  that  can  deliver  soft  organic  and  biologically  active  materials  onto  a  surface,  apply  a 
controlled  localized  force,  and  print  over  large  areas  with  micrometer  scale  resolution.  We  report 
how  arrays  composed  of  approximately  10,000  elastomeric  tips  mounted  on  the  piezoactuators 
of  an  atomic  force  microscope  (AFM)  can  be  used  to  apply  a  localized  force  that  accelerates  the 
Huisgen  1,3-dipolar  cycloaddition  reaction  between  azide-terminated  monolayers  and  alkyne 
inks  in  the  absence  of  the  Cul  catalyst.  Upon  bringing  the  tips  into  contact  with  the  surfaces,  a 
nanoreactor  is  formed  between  the  tips  and  the  surfaces  that  confines  the  reactants  under  an 
applied  force  (Figure  16).  Several  aspects  of  this  approach  are  ideal  for  patterning  applications 
and  for  understanding  how  force  and  interface  structure  contribute  to  the  rate  constant  of 
reactions,  namely  (1)  these  tip  arrays  uniformly  deposit  delicate  organic  and  biologically  active 
materials  onto  surfaces,  (2)  the  piezoelectric  tip  actuation  provides  precise  control  over  reaction 
time  and  the  applied  force  between  the  tips  and  the  substrates,  and  (3)  the  large  patterns  (cm^) 
created  by  the  massively  parallel  tip  arrays  facilitate  fluorescence  microscopy  and 
electrochemical  analysis  of  the  surface  reactions.  In  these  studies  we  show  that  the  rate 
constant  k  of  the  Huisgen  reaction  on  surfaces  is  sensitively  dependent  on  applied  force  and, 
surprisingly,  the  alkyl  chain  length  of  the  azide-terminated  monolayer. 

The  CV  studies  were  used  to  determine  the  activation  volume  \/*obs  (observed 
dependence  of  In  k  on  pressure)  by  plotting  In  k  versus  P.  An  activation  volume  of  -189  ±  41 
cm^  moF^  was  determined  for  the  Huisgen  cycloaddition  on  the  C11  1 1-azidoundecane-1 -thiol 
SAM.  Calculations  at  the  M06-2X/6-31G(d)  level  in  the  gas  phase  show  that  the  (3  -i-  2) 
cycloaddition  of  alkyl  azide  and  terminal  alkyne  has  a  V*  of  about  -13  cm^  mofV  This  value  is 
significantly  smaller  than  that  observed  in  the  experiment.  We  attribute  the  difference  between 
these  two  values  to  the  monolayer  compression  that  occurs  during  printing.  We  propose  that  the 
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Fig.  16  (a)  An  elastomeric  tip-array  coated  with  an 
ink  mixture  consisting  of  alkyne  and  PEG  (red) 
printing  onto  an  azide-terminated  Au  surface,  (b) 
Magnified  view  of  the  inked  tip  array  approaching 
the  surface,  (c)  The  tip  forms  a  nanoreactor  with 
the  surface  by  applying  a  force  that  accelerates 
the  Huisgen  cycloaddition  reaction,  (d)  After 
washing,  only  covalently  bound  molecules  remain 
on  the  surface. 


observed  activation  volume  V*obs  in  the  nanoreactor  is  the  sum  of  the  activation  volume  of  the 
cycloaddition  reaction  V^^c  and  the  change  in  the  volume  of  the  monolayer 

At/obs  =AV*+V^ 

To  test  this  hypothesis,  we  reasoned  that  the  height  of  the  monolayer  would  alter  V^^obs 
because  of  the  differences  in  the  change  of  A 14  during  compression.  Therefore,  2  was  printed 
on  SAMs  containing  alkane  chains  of  different  length  between  the  thiol  and  the  azide  with  a 
printing  time  of  5  min  and  at  pressures  of  2.72  -  7.95  MPa  (Figure  17).  Based  on  a  maximum 
packing  density  of  2.7  x  10^"^  molecules  per  cm^,  the  measured  \Aobs  decreases  gradually  from 
-189  cm^  mol"^  for  an  undecyl  chain  to  -39  cm^  mol”^  for  a  propyl  chain.  This  indicates  that, 
assuming  that  V*c  remains  constant  (-13  cm^  moF^),  the  volume  change  of  the  monolayer  Al/m 
becomes  smaller  and  smaller  with  the  decrease  of  the  height  of  the  monolayer.  A  simple  model 
can  semiquantitatively  account  for  the  change  in  A\4.  The  free  rotation  of  the  C11  azide 

molecule  on  Au  in  a  standing  up 
configuration36  sweeps  out  a  cylinder.  From 
the  density  of  2.7  x  10^"^  molecules/cm^  of 
the  Cl  1  chain  on  the  surface,  the  occupied 
area  of  the  single  molecule  is  calculated  to 
be  0.37  nm^.  The  height  of  the  Cl  1 
monolayer  is  about  2.0  nm,  so  the  volume 
of  the  cylinder  is  0.74  nm^  or  445  cm^  mofV 
Calculations  indicate  that,  when  two  anti- 
conformers  in  the  middle  of  the  alkane  chain 
are  changed  to  gauche  in  the  reaction,  the 
height  of  the  Cl  1  monolayer  decreases  to 
1 .2  nm.  This  will  result  in  a  Ah  (Figure  1 7)  of 
-0.8  nm,  leading  to  an  A\4  of  -0.30  nm^  or 
-178  cm^  mor\  According  to  this  model, 
the  computed  A\4  values  for  the  C8,  C6, 
C5,  and  C3  monolayers  are  -156,  -134, 
-111,  and  -78  cm^  mor\  respectively,  and 
these  data  reproduce  the  trend  observed  in 
the  experiments.  The  largest  discrepancy 
between  experimental  and  computational 
values  occurs  with  the  C3  monolayers, 
which  are  more  likely  to  lie  flat  on  the 
substrate  and  are  therefore  less 
compressible  along  the  force  vector  than 
the  closely  packed  monolayers  used  for 
computations. 

In  conclusion,  we  have  used 
massively  parallel  tip  arrays  to  form  nanoreactors  that  can  apply  force  to  accelerate  bond 
forming  reactions  with  negative  Vt  on  surfaces.  This  approach  was  used  to  show  that  the 
velocity  of  the  Huisgen  reaction  on  surfaces  is  sensitively  dependent  on  force  and  monolayer 
chain  length,  and  a  semiquantitative  model  was  developed  that  explains  the  previous 
contradictory  results  in  the  literature.  Force  may  have  an  important  role  in  many  interfacial 
processes,  and  the  approach  described  in  this  work  can  be  used  to  address  challenges  in 
materials  science,  biology,  and  nanotechnology.  The  results  of  these  efforts  are  described  in 
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Fig.  17  (a)  Plot  of  Ink  vs  Pat  different  alkyl  chain 
iengths,  whose  slope  provides  l^obs-  Error  vaiues 
were  determined  from  three  independent 
measurements,  (b)  The  rotation  of  the  Cl  1  azide 
molecule  sweeps  out  a  large  cylinder,  (c)  The 
monolayer  is  compressed  when  tips  are  brought 
into  contact  with  surface,  decreasing  the  cyiinder 
voiume  in  the  reaction  transition  state. 
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further  detail  in  an  article  recently  published  in  the  Journal  of  the  American  Chemical  Society  {J. 
Am.  Chem.  Soc.  2014,  136,  10553). 
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Fig.  19  While  self-propelled  rods  are 
trapped,  they  move  at  speeds 
comparable  to  when  they  are  far 
away  from  any  sphere.  Symbols  and 
error  bars  respectively  show  the 
temporal  mean  and  one  standard 
deviation  in  speed  of  each  rod.  Solid 
and  dashed  lines  have  slopes  0.9  ± 
0.1. 


Fig.  18  Sample  trajectories  of  self- 
propelled  rods  moving  at  speeds  in  the 
order  of  20  pm/s  near  passive 
spheres,  (a)  A  rod  encounters  a 
sphere,  orbits  around  it  five  times,  and 
moves  away,  (b)  A  succession  of 
orbits  with  the  same  sense  of  rotation, 
(c)  The  sense  of  rotation  switches,  (d) 
A  slalom  through  a  series  of  spheres. 


C.5.  Catalytically  Propelled  Nanorod  Dynamics 

A  major  aim  of  the  proposed  work  was  to  couple  inorganic 
nanostructures  with  the  artificial  carbohydrate  receptors  to 
develop  synthetic  systems  with  dynamics  and 
computational  complexity  similar  to  biological  networks.  In 
Y1,  we  reported  on  a  new  set  of  equations  that  describe 
quantitatively  the  motion  of  active  catalytic  nanorods  in 
H2O2  solution  {Physics  Review  Letters,  2013,  110, 

038301).  These  equations  provide  a  predictive  foundation 
for  understanding  the  motion  of  real  rods  possessing 
defects  that  are  the  inevitable  result  of  any 

nanomanufacturing 
process,  and 

therefore  accurately 
predicted  the 

dispersion  of  activated 
nanorods  by 

accounting  for  flipping 
that  occurred  as  a 
result  of  these 
defects.  In  Y2  we 
extended  this  work  by 

investigating  the  interactions  of  these  active  AuPt  nanorods 
with  spherical  microparticles  to  determine  how  barriers 
alter  the  motion  of  these  microswimmers,  thereby  providing 
a  more  complete  description  of  their  motion. 

We  show  that,  in  a  suspension  of  passive  spheres 
on  a  horizontal  surface,  self-propelled  Au-Pt  rods  are 
temporarily  trapped  and  orbit  closely  around  the  spheres 
(Figure  18).  These  observations  are  consistent  with  our 
model  based  on  lubrication  theory,  which  predicts  that 
swimmers  are  attracted  to  obstacles  over  short  distances. 
The  only  key  ingredient  in  the  model  is  a  prescribed  slip  on  the  swimming  surface,  a  suitable 
assumption  for  modeling  the  motion  of  a  variety  of  phoretic  swimmers  and  motile  ciliates.  Thus, 
the  attraction  should  have  broad  implications,  ranging  from  the  applications  of  swimmers  in 
synthetic  devices  to  the  study  of  cell  adhesion  and  feeding.  We  fabricate  Au-Pt  rods  with  length 
L  =  2  ±  0.2  pm  and  diameter  0.39  ±  0.04  pm  following  the  method  of  electrochemical  deposition 
in  anodic  aluminum  oxide  membranes.  The  rods  and  passive  spheres  with  typical  diameter  5.5 
±  0.5  pm  are  suspended  in  a  H2O2  solution  with  typical  concentration  15%.  The  system  here  is 
effectively  two  dimensional  because  the  rods  and  spheres  remain  close  to  the  microscope  slide. 
The  positions  of  the  rods  and  spheres  are  tracked  using  optical  microscopy  (Nikon  Eclipse  80i, 
40X),  digital  camera  (Lumenera  Infinity  1-3),  and  software  (Imaged  and  Matlab).  In  the  absence 
of  passive  spheres,  the  rods  turn,  flip,  and  disperse  over  time  as  reported  previously.  The 
presence  of  passive  spheres  significantly  alters  the  trajectories  of  self-propelled  rods  (Figure 
18).  When  a  self-propelled  rod  encounters  a  sphere,  the  rod  typically  orbits  around  it.  The  rod 
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Fig.  20  Image  sequence  of  a  solid  sphere  which  hardly 
rotates  while  a  self-propelled  rod  orbits  around  it.  The 
sphere  undergoes  minor  translational  diffusion  confirming 
that  it  is  not  fixed  to  the  substrate.  Red  arrows  point  to  a 
tracer  on  the  sphere.  Blue  dots  and  arrows  show  the 
location  and  direction  of  motion  of  the  rod. 


can  move  around  a  succession  of 
spheres  and  either  continue  to  turn  in  the 
same  sense  or  switch  handedness  of 
circular  motion.  The  trajectory  may  deflect 
only  slightly  when  the  rod  is  close  to  a 
sphere,  but  the  direction  is  always 
towards  the  sphere. 

This  phenomenon  of  self- 
propelled  rods  orbiting  around  spheres 
occurs  over  a  wide  range  of  materials  and 
sizes.  We  observe  rods  ranging  in  length 
from  1  to  4  pm  and  spherical  beads  made  of  glass,  polymerized  3-(trimethoxysilyl)propyl 
methacrylate  (TPM),  poly(methyl  methacrylate)  (PMMA),  and  polystyrene,  ranging  in  diameter 
from  1  to  125  pm.  Despite  differences  in  chemical  and  physical  properties  on  their  surface,  the 
spheres  can  trap  self-propelled  rods  in  all  cases,  except  for  1  pm  polystyrene  beads.  This 
exception  is  attributed  to  the  large  sedimentation  height  of  the  polystyrene  beads  which  do  not 
remain  close  to  the  rods.  The  self-propelled  rods  almost  maintain  their  typical  speed  while  they 
temporarily  orbit  around  spheres  (Figure  19).  This  shows  that  the  rods  hardly  experience  any 
additional  drag  despite  moving  close  to  the  solid  spheres.  While  a  rod  is  trapped,  it  remains  in  a 
narrow  region  between  the  sphere  and  the  surface  (Figure  20),  whose  height  is  estimated  to  be 
-550  nm  after  accounting  for  minor  optical 
aberrations.  This  height  is  not  much  greater 
than  the  rod  diameter,  meaning  that  the  rods 
orbit  closely  around  spheres. 

Do  spheres  rotate  while  a  rod  orbits 
around  them?  To  answer  this  question,  we 
use  inert  PMMA  spheres  with  an  embedded 
piece  of  hematite  acting  as  a  black  marker. 

The  marker  remains  almost  stationary  while 
a  rod  orbits  around  a  sphere  (Figure  21).  The 
sphere  is  almost  always  detached  from  the 
moving  rod,  suggesting  that  any  collision 
occurs  only  briefly.  This  is  the  first 
experimental  evidence  supporting  the 
previous  hypothesis  that  the  rods  generate  a 
local  flow  along  their  surface,  as  elaborated 
further  in  our  theory  below.  There  are  several 
possible  mechanisms  of  trapping.  The  rods 
are  slightly  curved  and  tend  to  move  in 
curved  paths  with  typical  curvature  k~  0.12 
pm'^  so  they  could  slide  along  at  and  curved 
surfaces  with  curvature  much  less  than  k. 

However,  we  observe  rods  orbiting  around 
glass  spheres  with  diameter  as  small  as  1 
pm,  which  has  curvature  much  greater  than 
K.  Long-range  hydrodynamic  attraction 
towards  a  surface  is  possible  for  swimmers  that  are  represented  by  a  force  dipole  with  positive 
strength  (pusher),  though  our  rods  experience  only  short-range  attraction  towards  spheres  as 
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Fig.  21  (a)  Sketch  of  a  rod  orbiting  around  a  sphere  on 
a  horizontal  surface,  (b)  Theoretical  configuration  of  a 
swimming  body  near  a  rigid  wall  in  the  reference  frame 
of  the  swimmer.  The  wall  approximates  the  local 
surface  of  the  sphere  and  translates  to  the  right.  (c,d) 
Predictions  for  the  swimming  speed  U,  distance  hO 
away  from  the  wall,  and  rescaled  slope  m  of  a  swimmer 
with  constant  slip  velocity  prescribed  on  90%  of  the 
front  or  back  of  the  body.  Dashed  lines  show  the  early- 
time  regime. 


confirmed  below.  The  attraction  could  be  a  consequence  of  short-range  hydrodynamic 
interactions  much  like  between  bacteria  and  walls.  The  short-range  interactions  between  rods 
and  spheres  are  studied  further  by  formulating  a  two-dimensional  model  of  a  rod  moving 
horizontally  near  a  wall.  While  our  experimental  system  is  three  dimensional  with  curved 
boundaries,  our  simplified  model  approximates  the  local  surface  of  the  spherical  obstacle  by  a 
flat  immobile  wall.  Thermal  fluctuations  are  neglected  for  simplicity.  We  adopt  a  Cartesian 
coordinate  system  (x,  y)  in  the  reference  frame  that  translates  with  the  swimmer  along  the  wall, 
where  x  is  scaled  by  a  half  of  the  length  of  the  swimming  body  L=2  and  y  is  scaled  by  a  typical 
distance  H  between  the  swimmer  and  the  wall  (Fig.  16b).  In  this  moving  reference  frame,  the 
wall  is  represented  by  y  =  0  and  translates  to  the  right  with  speed  U,  where  all  speeds  are 
scaled  by  a  typical  slip  velocity  V  on  the  swimmer.  The  wall-facing  surface  of  the  swimmer  is 
represented  by  y  =  h  with  h  =  ho(1  -i-  mx),  where  ho  is  a  re-scaled  distance  from  the  wall  to  the 
mid-point  of  the  swimming  surface  and  m  is  a  re-scaled  slope  of  the  swimmer  relative  to  the 
wall.  Note  that  m  =  1  corresponds  to  the  front  end  of  the  swimmer  hitting  the  wall.  The  aim  is  to 

predict  how  U,  ho,  and  m  evolve 
over  time,  given  a  prescribed  slip 
speed  us  on  the  surface  of  the 
swimmer. 

To  verify  experimentally  that 
only  short-range  interactions  occur 
between  self-propelled  rods  and 
passive  spheres,  we  consider  the 
effective  density  of  a  suspension  of 
rods  as  a  function  of  the  distance 
away  from  the  lowest  point  of  a 
sphere  (Figure  22).  This  is 
measured  by  sampling  the  distance 
from  a  specified  rod  to  the  nearest 
sphere  over  time,  excluding  cases 
when  the  rod  has  multiple  spheres 
within  three  radii  of  the  sphere.  For 
passive  rods  in  water,  the  density  is 
almost  uniform  everywhere  except 
below  the  sphere,  where  the  density 
decreases  to  zero.  For  active  rods 
in  H2O2,  the  density  is  almost 
uniform  everywhere  except  below  the  sphere  again,  where  the  rods  are  now  highly 
concentrated.  This  supports  our  theory  that  moving  rods  remain  trapped  primarily  because  of 
short-range  hydrodynamic  effects. 

To  explore  why  the  rods  eventually  move  away  from  the  spheres,  we  observe  that  the 
trapping  times  are  typically  less  than  5  seconds  and  that  their  distribution  shows  an  exponential 
decay  with  rate  A=1.5  s'^  (Figure  22).  In  rare  instances,  the  rods  remain  trapped  for  longer  times 
in  the  order  of  minutes.  Flighly  curved  rods  that  move  in  small  circles  tend  to  either  remain 
trapped  for  relatively  long  times  or  hardly  remain  trapped,  depending  on  whether  they  approach 
a  sphere  with  the  trajectory  curving  towards  or  away  from  the  sphere,  respectively.  Perhaps  the 
rods  are  curved  towards  the  sphere  while  they  are  trapped  until  fluctuations  flip  them  over  and 
enable  them  to  escape,  suggesting  that  rod  curvature  and  flipping  affect  the  movement  of  rods 
near  spheres.  Thermal  fluctuations  are  most  likely  to  be  driving  the  rods  away  from  remaining 
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Fig.  22  (a)  Experimental  data  of  the  effective  density  of  a 
suspension  of  rods,  scaled  by  the  expectation  if  the  rods  were 
distributed  uniformly.  This  is  plotted  against  the  distance  from  the 
lowest  point  of  a  sphere  in  units  of  the  sphere  radius,  based  on 
over  4000  measurements  of  the  distance  over  time.  Below  the 
sphere,  the  density  of  passive  rods  falls  while  the  density  of  self- 
propelled  rods  peaks.  The  excluded  region  is  where  the  rods  are 
physically  unable  to  enter,  (b)  Histogram  of  trapping  times  of  self- 
propelled  rods  orbiting  around  spheres.  Inset  shows  the 
probability  that  the  rods  remain  trapped  beyond  a  specified  time. 
Red  curves  are  exponential  distributions  with  rate  parameter  A  = 

1.5  s'\ 


trapped,  since  the  time  scale  associated  with  trapping  is  comparable  to  the  time  scales  of 
rotational  diffusion  and  duration  between  stochastic  flips,  which  are  both  in  the  order  of  one 
second.  While  our  deterministic  model  predicts  that  swimmers  move  closer  to  a  confining 
surface,  fluctuations  may  cause  the  swimmers  to  move  away  or  remain  at  a  characteristic 
distance  from  the  surface.  This  may  explain  why  self-propelled  rods  appear  to  move  closer  to 
the  substrate  than  the  expected  sedimentation  height  of  purely  diffusive  rods.  The  exact 
distance  between  swimmers  and  confining  boundaries  could  be  predicted  in  the  future  by 
extending  our  model  to  incorporate  the  effects  of  fluctuations,  gravity,  and  the  curvature  of  the 
swimmers  and  boundaries. 

The  short-range  attraction  towards  walls  has  broad  implications  on  potential 
applications  of  synthetic  swimmers  as  well  as  addressing  the  behavior  of  motile  cells.  Many 
swimmers  including  motile  ciliates  and  synthetic  phoretic  particles  could  be  interpreted  to  have 
an  effective  slip  on  a  majority  of  their  surface.  Our  model  suggests  that  these  swimmers 
accumulate  along  boundaries,  contrary  to  passive  particles  which  have  difficulty  diffusing  into 
spatially  constricted  regions.  This  may  offer  insight  into  the  possible  hydrodynamic  effects  in  cell 
adhesion  and  feeding,  as  well  as  suggest  strategies  to  capture  and  filter  groups  of  swimmers 
using  synthetic  devices.  Our  model  also  implies  that,  while  the  swimmers  may  not  mix  fluids 
over  long  distances,  payloads  attached  to  the  swimmers  may  exert  more  drag  and  perhaps 
develop  large-scale  flow  patterns  as  observed  in  bacterial  baths.  Short-range  interactions 
between  swimmers  and  obstacles  importantly  affect  their  movements  and  must  be  incorporated 
in  any  accurate  model  of  self-propelled  colloidal  suspensions,  and  will  be  particularly  important 
when  using  these  systems  to  understand  how  carbohydrate  recognition  can  influence  the 
assembly  of  dynamic  nanostructures.  The  results  of  these  efforts  are  described  in  further  detail 
in  an  article  recently  published  in  Soft  Matter  {Soft  Matter  20^^,  10,  1784). 


C.6.  Summary  of  Accomplishments 

In  summary,  this  award  from  the  AFOSR  has  directly  supported  several  major  research 
accomplishments  that  will  lead  to  the  creation  of  carbohydrate-based  nanotechnology. 
Specifically,  we  have  created  a  synthetic  lectin  that  employs  cooperativity  to  obtain  selectivity 
for  mannose  and  characterized  the  thermodynamics  of  binding  and  the  structure  of  this  receptor 
when  bound  to  several  monosaccharides.  This  work  is  a  major  breakthrough  because,  white 
cooperativity  is  ubiquitous  in  the  binding  of  carbohydrates  in  nature,  our  work  is  the  first 
example  that  demonstrates  cooperativity  in  a  synthetic  carbohydrate  receptor.  Another  major 
accomplishment  is  that  we  have  developed  new  tools  and  techniques  to  create  glycan  arrays, 
including  a  photochemical  approach  that  leads  to  glycopolymers  that  are  orders  of  magnitude 
more  sensitive  towards  lectins  than  glycan  arrays  prepared  by  other  methods.  In  the  course  of 
these  studies,  we  have  also  increased  the  understanding  of  reactivity  of  surfaces  in  several 
important  ways,  including  developing  the  first  method  to  covalently  pattern  graphene  through  a 
force-accelerated  Diels-Alder  reaction  and  also  developed  a  new  quantitative  understanding  of 
how  force  affects  reactivity  on  self-assembled  monolayers.  Finally,  we  have  developed  a 
quantitative  description  for  the  motion  of  catalytic  Au-Pt  nanoswimmers  that  explain  their 
diffusion  and  interaction  with  barriers.  Taken  together,  these  results  demonstrate  the  potential  of 
carbohydrates  as  directing  elements  in  advanced  nanotechnologies  and  provide  a  strong 
foundation  for  developing  this  transformational  research  area. 
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